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ABSTRACT 

Context. Most stars are born in clusters, thus the protoplanetary discs surrounding the newly formed stars might be influenced by this 
environment. Isolated star-disc encounters have previously been studied, and it was shown that very close encounters are necessary to 
completely destroy discs. However, relatively distant encounters are still able to change the disc size considerably. 

Aims. We quantify the importance of disc-size reduction that is due to stellar encounters in an entire stellar population. 

Methods. We modelled young, massive clusters of different densities using the code Nbody6 to determine the statistics of stellar 
encounter parameters. In a second step, we used these parameters to investigate the effect of the environments on the disc size. For 
this purpose, we performed a numerical experiment with an artificial initial disc size of 10 5 AU. 

Results. We quantify to which degree the disc size is more sensitive to the cluster environment than to the disc mass or frequency. We 
show that in all investigated clusters a large portion of discs is significantly reduced in size. After 5 Myr, the fraction of discs smaller 
than 1000 AU in ONC-like clusters with an average number density of J> clmKr ~ 60 pc -3 , the fraction of discs smaller than 1 000 AU 
is 65%, while discs smaller than 100 AU make up 15%. These fractions increase to 84% and 39% for discs in denser clusters like IC 
348 (fl c i us t er ~ 500 pc -3 ). Even in clusters with a density four times lower than in the ONC (p dusler ~ 15 pc -3 ), about 43% of all discs 
are reduced to sizes below 1 000 AU and roughly 9% to sizes below 100 AU. 

Conclusions. For any disc in the ONC that initially was larger than 1 000 AU, the probability to be truncated to smaller disc sizes as a 
result of stellar encounters is quite high. Thus, among other effects, encounters are important in shaping discs and potentially forming 
planetary systems in stellar clusters. 

Key words. Methods: numerical - Protoplanetary disks - Stars: kinematics and dynamics - Galaxies: star clusters: general 


1. Introduction 


Table 1: Observed disc radii in different cluster environments. 


Most, if not all stars are initially surrounded by a protoplanetary 
disc from which a planetary system might eventually form. The 
observed sizes (= radii) of these discs span from a few 10 AU 
to several 100 AU, see Table Q](Vicen fe~& Alves||2005 Eisner 


et al. 2008[ |McC aughrean & Q’de ll|1996 Andrews & Williams 


2007 Andrews et al.|2009||2010| . However, the majority of disc 
sizes has been found to be in the range of 100 - 200 AU. These 
discs disperse rapidly as a result of internal processes such as 
torques (e.g. Shu et al. 19871, turbulent effects (|Klahr 


& Bodenheimer 2003], and magnetic fields (Balbus & Hawley 


2002). Observations of disc fractions in clusters show that the 
discs dissipate through internal processes on average at an age of 
< 2.5 Myr. It is currently unclear whether this is valid for young 
stars in general ( jPfalzner et al.||2014] and how certain absolute 
age estimates in clusters are ( Bell et al.|2013] . 

There are additional environmental processes in a cluster 
that, depending on its stellar density, might lead to disc-mass 
loss, reduction in physical size, or even complete destruction. 


These processes are photoevaporation (see e.g. Johnstone et al. 


1998 )|Storzer & Hollenbach|1999||Scally & Clarke|2002j |Mat 


suyama et al.||2003| |Johnstone et al.||2004[ Adams et al.||2006| 


Alexander et al. 20061 Ercolano et al.||20081 IGorti & Hollen^] 

bach|2009 Drake et al.|2009] and stellar fly-bys or gravitational 


T993||Scally & Clarke|2002l 


stripping (see e.g. Clarke & Pringle 


Pfalz ner|20041 1 Adams et al.||2006T Olczak et al.|2006||Pfalznet1 


Cluster 

Disc radii [AU] 

N sources 

References 

ONC 

50 - 200" 

10 

(1) 

ONC 

80-217 

39 

(2) 

ONC 

27 - 506“ 

6 

(3) 

Ophiuchus 

25 - 700 

11* 

(4) 

Ophiuchus 

14-198 

16* 

(5), (6) 

Ophiuchus 

165 - 190 

2 

(7) 

Taurus 

25 - 600 

12 

(4) 

Taurus 

< 50 - 100 

4 

(8) 


Notes. Observed clusters (Col. 1) and disc-radius ranges (Col. 2). The 
number of sources V sources and the references are given in Cols. 3 and 4. 

Values calculated from given disc diameters. <A) The sources in Cols. 
(4), (5), and (6) overlap. 


References. (1) Vicente & Alves (2005); (2)|Eisner et al. 

(2008); (3) 

McCaughrean & O’dell (1996); (4) Andrews & Williams 

(2007); (5) 

Andrews et al. (2009); (6) Andrews et al. (2010); (7) Brinch & Jpr- 

gensen (2013); (8) Harsono et al. (2014). 


& Qlczak||2007[ |Qlczak et al ||2010| [Craig & Krumholz||20131 


Steinhausen & Pfalzner|2014||Rosotti et al.|2014[). 

From observations, the sizes of protoplanetary discs are 
usually determined by fitting the spectral energy distributions 


(SEDs) to disc models using a truncated power-law (Andrews 
|& Williams||2007[). If resolved images of the discs can be ob- 
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tained, the disc size is identified as the radius at which the lu¬ 


minosity drops below a certain value (Vicente & Alves 2005 
O’dell 1998). However, it is much more difficult to observation- 
ally constrain the disc size than the disc frequency or disc mass. 
Therefore, the disc sizes are only determined for a very limited 
number of stars. As a consequence, it is still a subject of debate 
whether the disc size is a function of the mass of the central star 
(compare e.g. Hillenbrand et al. 1998} Vicente & Alves| 2005 


Eisner et al.|2008 1 . The disc frequencies in clusters of similar age 


(~ 1 Myr) reveal that the cluster environment influences the pro¬ 
toplanetary discs. The percentage of stars surrounded by discs 
decreases with increasing cluster peak density, dropping from 
85% in NGC 2024 and 80% in the Orion nebula cluster (ONC) 
down to 40% in NGC 3603 (see also Table |2}. Because the clus¬ 
ter ages are similar, the higher cluster density leads to a higher 
degree of disc destruction and thus to a lower disc frequency. 

A number of theoretical investigations studying the influence 
of individual stellar encounters on protoplanetary discs have 
been performed in the past. For the case of an equal-mass en¬ 
counter, Clarke & Pringle ( 1993} were the first to study the fate 
of disc matter after an interaction. They found that a coplanar, 
prograde encounter at a relative periastron distance of 1.25r pcn 
material is removed from the disc down to 0.5r pL , n . But mate¬ 
rial is not only removed, it is also re-distributed within the disc, 


changing the (mass) surface density profile of the disc (Hall 
[1997). As the disc masses of individual stars are observation- 
ally better constrained than disc sizes, there has been a consid¬ 
erable amount of theoretical research on the influence of stellar 


encounters on the disc mass (e.g. Olczak et al. 


2010 


& Pfalzner 2014), but fewer studies of the disc size (e.g. |Hall 


Steinhausen 


et al.|1996| Kobayashi & Ida|2001[ Breslau et al.|2014 ). Never¬ 


theless, the disc size is an important parameter because it defines 
the physical extent of the potentially forming planetary system. 

The disc size after an encounter event was investigated for 
example by |Kobayashi & Ida| ( [2001 j ). They found that the disc 
size after an equal-mass encounter is one-third of the periastron 
distance. Adams ( 2010| ) found a slightly different disc size of b/3 
where b is the impact parameter, which is not necessarily equal 
to the periastron distance. Two other estimates for disc-size up¬ 
per limits were introduced by de Juan Ovelar et al. ( j2012[ ). They 
derived the disc size assuming that all material outside the La- 
grangian point between the two stars is stripped from the disc. 
Furthermore, they suggested to convert the formula for disc- 
mass loss for parabolic, coplanar encounters reported by Qlczak| 


et al. (2006) directly into an upper limit for disc truncation, as¬ 
suming a mass-density distribution within the disc of Z = r~ l . 
The few previous studies that considered disc sizes in stellar 
cluster environments (Adams et al. 2006, Malmberg et al. 2011; 


Pfalzner 


2013|) generalised, for lack of a universal description 


of disc sizes, the results obtained for equal-mass encounter by 
Kobayashi & Ida (2001)) to non-equal mass encounters. 


Recently, Breslau et ah] ( |2014| ) showed that a generalisation 
of this one-third criterion postulated by Kobayashi & Ida ( 2001| ) 
is not advisable and that the outcome is very sensitive to the mass 
ratio of the two stars. They performed simulations of isolated 
star-disc encounters for a wide range of periastron distances and 
mass ratios, choosing the highest contrast in surface density to 
determine the disc sizes after the encounters. This is similar to 
the criterion used by observations, and the results by |Breslau| 


et al. (2014) can easily be applied to cluster simulations. They 


showed that encounters with large periastron distances and/or 
low mass ratios do not remove material from the disc but shift it 
inwards, reducing the disc size but not the mass. For disc-mass 
loss closer encounters and/or larger mass-ratios are needed. This 
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Table 2: Observed cluster properties. 


Cluster 

p[ pc~ 3 ] 

f&isc 

References 

NGC 2024 

* 10" 

85% 

(1), (2) 

ONC 

* 10 4 

80% 

(3), (4) 

NGC 3603 

* 10 5 

40% 

(5), (6) 


Notes. Col. 1 gives the cluster name, p is the cluster density (Col. 2), 
and /disc the disc fraction (Col. 3). The references for the densities and 
disc fractions are given in Col. 4. 


References. (1) Lada et al. 

(1991); (2) Haisch et al. 

2000); (3) Me- 

Caughrean & Stauffer 

(1994 

i; (4)]Lada etal. (2000); (5) 

Harayama et al. 

(2008 ); (6) Stolte et a 

. (2004). 


finding demonstrated that the upper disc-size limits of |de Juan| 
Ovelar et al. ( 2012| >, which connect disc-size change directly to 


disc-mass loss, are thus too high by up to a factor of two. 

Recently, [Rosotti et al. ( 2014[ ) investigated the effect of en¬ 
counters on the disc-size distribution in an entire stellar clus¬ 
ter. They performed combined smoothed particle hydrodynam¬ 
ics (SPH)/N-body simulations of a low-mass cluster (100 stars) 
allowing hydrodynamic al disc spreading. They found that while 
the masses of the discs are not much affected, the disc sizes are 
reduced by encounters within the ~ 0.5 Myr of cluster evolution. 

Here, we wish to analyse numerically how encounters influ¬ 
ence the disc size in clusters using a realistic initial mass function 
(IMF) and the therefore broad spectrum of encounter mass-ratios 
found in real clusters. We model clusters of different density to 
quantify the resulting disc-size dependence. Most importantly, in 
contrast to previous work, the effects of the encounters are mod¬ 
elled taking into account that the final disc size is very sensitive 
to the mass ratio of the two interacting stars. Although exter¬ 
nal photoevaporation might be as important as stellar encounters 
when describing disc sizes in clusters, we do not consider it in 
this work. 

The simulations of the clusters with different densities as 
well as the disc-size definition are introduced in Sect. [2] The 
choice of initial disc sizes and the classification of encounters 
used throughout this work are also discussed there. The most 
important results are presented in Sect. [3 the assumptions ap¬ 
plied and their influence on these results are detailed in Sect. |4] 
Our findings are briefly summarised in Sect. [5] 


2. Method 


We performed numerical simulations to investigate the influence 
of star-disc encounters on the disc size within stellar clusters 
of different densities. This was achieved by performing Nbody6 
simulations ( Aarseth|1973[ Spurzem||1999| 


ilar to those of|Qlczak et al. (|2008| ), |Qlczak et al.| (2010 1 , and 


Aarseth 


2003 i sim 


|Steinhausen & Pfalzner ( 2014 ). In those simulations, the stel¬ 
lar encounters in the clusters were tracked and their effect on 
protoplanetary discs was analysed in the diagnostics using a pre¬ 
viously performed parameter study. In contrast to the previous 
work, we investigate the disc size rather than the disc mass or 
frequency in clusters. The disc sizes were obtained using the re¬ 
sults from lBreslau et al.[ ( |2014| >. 


2.1. Cluster simulations 

The dynamics of the stellar clusters was modelled by Nbody6 
simulations taking the cluster density profile of the ONC as a 
prototype for a massive cluster in the solar neighbourhood. The 
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ONC is the best observed massive cluster, rendering it possible 
to compare results of theoretical studies to observational data. 
It contains about 4000 stars, is ~ 1 Myr old, has a half-mass 
radius of roughly 1 pc (Hillenbrand & Hartman n) 1998|>, and ex¬ 
tends to ~ 2.5 pc ([Hillenbrand & Hartmann||1998| |Hillenbrand 


& Carpenter |2000 1 . Lada et al.| ( |2000j l found that in the Trapez 

ium, that is, in the inner 0.3 pc of the ONC, about 80 - 85% of 
the stars are surrounded by a protoplanetary disc, 60% of which 
having radii smaller than 50 AU ( Vicente & Alves|2005| ). 

|Scally et ~aI7| ( |2005) > found the density distribution in the ONC 
to differ from the isothermal approximation po CT >' 2 in the inner 
0.2 pc of the cluster. The surface density there is much flatter 
To obtain this flat profile after 1 Myr of simula- 


with X oc 


„- 0.5 


tion time, which is the assumed age of the ONC, Hillenbrand 
( |1997| > and |Qlczak et al.| ( |2010| >, for instance, used the three- 
dimensional initial number density distribution 


PoO) = 


Po (r/flo. 2 ) 23 , r 6 (0, P 0 . 2 ] 
' Po ( r /R().2) - '° I r 6 (^0.2-^] 


0, r e (R , oo], 


(1) 


with the radius Rq.i = 0.2 pc, the cluster radius R = 2.5 pc, 
and po = 3.1 x 10 3 stars pc~ 3 . After 1 Myr of cluster evolution 
this distribution represents today’s number density distribution 
of the ONC. The stellar positions, masses, and velocities were 
randomly attributed to the stars. All stars were assumed to be 
initially single and without any stellar evolution throughout the 
simulation. This is standard procedure to reduce the computa¬ 
tional cost. Furthermore, the clusters were modelled without gas. 
The assumptions made and their influence on the outcome is dis¬ 
cussed in Sect.|4]in detail. 

The stellar masses were sampled from the IMF by [Kroupa| 

( 120021 


£(m) = 


m 1 ,0.08 < m/ M Q < 0.50 
, 0.50 <m/ M 0 < 1.00 
~ 2 ' 3 ,1.00 < m/ M 0 < oo, 


- 2.3 


( 2 ) 


with an upper mass limit of 150 M 0 and randomly attributed 
to the stars in the cluster. We neglected the effect of initial mass 


segregation as proposed for the ONC for example by Bonnell 
|& Davies] ( |1998| > and |Qlczak et al.| ( |2011) >. Using an IMF is 
important, as the effect of unequal-mass encounters onto discs 
is stronger than for equal-mass encounters (jOlczak et al.|2006[ 
|Moeckel & Bally|2007) . 

The clusters were set up to be initially in virial equilibrium: 


Q = = 0.5, where P hm is the half-mass radius, G the grav¬ 

itational constant, M the total cluster mass, and cr the velocity 
dispersion. The latter was sampled adopting a Maxwellian dis¬ 
tribution. 

As we are interested in the effect of clusters of different den¬ 
sities on the disc-size distribution, we set up density-scaled ver¬ 
sions of the ONC by keeping the radius fixed at 2.5 pc and vary¬ 
ing the number of stars within the cluster. The resulting clusters 
with 1000, 2 000, 8 000, 16 000, and 32 000 stars have average 
number densities of roughly 1/4, 1/2, 2, 4, and 8 times the den¬ 
sity of the ONC, respectively (see Tab lector the exact number 
and mass densities). 

For each of these clusters, a whole campaign of similar sim¬ 
ulations with different random seeds was performed to minimise 
the effect of initial stellar properties on the encounter history and 
thus the protoplanetary discs. The numbers of simulations A' s i rn 


Table 3: Cluster model set-up. 


Model 

N stars 

N sim 

Pcore 

[10 3 pc- 3 ] 

Pcluster 

[ pc~ 3 ] 

DO 

1000 

392 

1.3 

15.3 

D1 

2 000 

260 

2.1 

30.6 

D2 

4000 

264 

5.3 

61.1 

D3 

8 000 

78 

10.5 

122.2 

D4 

16 000 

21 

21.1 

244.5 

D5 

32 000 

14 

42.1 

489.2 


Notes. Model designation (Col. 1), initial number of stars in the cluster 
N saxs (Col. 2), number of simulations in the campaign N s i m (Col. 3), 
average number density in the cluster core p core (Col. 4), and average 
number density in the whole cluster p duster (Col. 5). 

(a) The cluster core radius was assumed to be the radius of the Trapezium 
(f^core = 0.3 pc). <A) The cluster radius is 2.5 pc. 

References. From |Qlczak et al,| i ]20 1 0[ i, [Stein hausen & Pfalzner) (2014) . 


per model are given in Table [3] This approach renders it possible 
to obtain good statistics of encounters and the resulting changes 
in the disc size for environments of different densities. 

In each simulation the encounter parameters for star-star in¬ 
teractions were recorded. In a diagnostic step, the disc sizes after 
each encounter were then calculated with these parameters, as 
described below in Sect. 12.21 


2.2. Encounter simulations 


In the diagnostic step the influence of each tracked encounter 
on the size of the protoplanetary discs was determined. To do 
this, we used the results of the parameter study by Breslau 
et al. ( 2014] ). Since their study was limited to prograde copla- 
nar parabolic (eccentricity= 1) encounters, which are the most 
destructive, the resulting disc sizes should be regarded as lower 


limits compared to hyperbolic and/or inclined encounters. Bres¬ 
lau et al. ( 2014j ) covered the parameter space for of mass ra¬ 
tios and periastron distances typical for star-disc encounters in 
cluster environments: m \2 — 0.3 - 500 and r p = 0.4 - 34, where 
m 12 = im/mi is the ratio between the perturber mass (m 2 ) and 
the mass of the disc-hosting star (mi) and r p = r pei i/Venous is 
the periastron distance r per i normalised to the disc size before 
the encounter r prev j ous . 

In a diagnostic step, we analysed the cluster simulations de¬ 
scribed above, using their fit formula 


Tlisc — 


0.28 ■ V peri 
/"previous 5 


0.: 

12 


if Liisc ^ ^"previous 
if Tji.sc — ^"previous i 


(3) 


with r per j being the periastron distance in AU, to calculate 
the disc size after each encounter in these clusters. As Eq.[3] 
shows, the resulting disc size is independent of the initial disc 
size as long as it is smaller than this initial value. It is basically 
the strongest encounter that determines the final disc size. At 
the beginning of the simulations, all disc sizes in the clusters 
were set to a universal value r; n ; t (for details see Sect. 13.21) . Most 
stars in a cluster undergo multiple encounters during the first few 
Myr. To take those multiple interactions of one and the same star 
with (different) other cluster members into account, the value of 
^previous for the star’s i th encounter was set to the disc size result¬ 
ing from the former (i - l)'* encounter (rp revious = r^,). This is 
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relative periastron distance encounter mass ratio 


Fig. 1: a) Relative final disc size (asterisks) and disc mass (open squares) after an equal-mass encounter event (mn = 1) as a 
function of the relative periastron distance r p = fp er i/rp r e V ious, with r per i periastron distance [AU] and r prev j ous disc size before the 
encounter [AU], b) Parameter pairs (encounter mass ratio mn and relative periastron distance r p ) leading to 5% reduction in disc- 
size (asterisks) and disc mass (open squares). The light and dark grey areas depict the parameter space in which the disc size and 
mass are reduced. 


an approximation as the discs might be “hardened 11 by stellar fly¬ 
bys or have enough time to spread again between two encounters 
due to viscosity (see Sect. |4]». The disc-size change due to multi¬ 
ple encounters needs to be investigated in more detail, taking all 
these effects into account. 

Additionally, we assumed that every star in the clusters was 
initially surrounded by a thin disc of a given size. Although 
this results in disc-disc encounters, the stellar interactions were 
treated as star-disc encounters using Eq. [3] The resulting disc 
sizes are not expected to differ much from the sizes after disc- 
disc encounters. For a detailed discussion of the approximations 
above, see Sect.[4]and |Breslau et al.[ ( |20T4| l. 

Before we present the results, we here clarify some defini¬ 
tions. The code Nbody6 that we used to simulate the clusters was 
extended to track the stellar encounters and store their properties. 
For every time step the star with the strongest gravitational in¬ 
fluence on a primary star was identified and its orbit calculated, 
for additional restrictions and details see |Qlczak et al.| ( |2Q12| >. 
Most such ’’encounters 11 are very distant and do not influence the 
sizes of discs at all. On the other hand, in extreme cases, stel¬ 
lar encounters cannot only alter the disc size, but destroy them 
completely. Here, the notion “disc-size-reducing encounter 11 is 
used for any encounter that (a) reduces the disc size by at least 
5% relative to its previous disc size and (b) does not destroy 
the disc (rd; sc > 10 AU). Consequences of this disc-size reduc¬ 
tion limit are discussed in Sect. [4] By contrast, the first en¬ 
counter that reduces the disc size by at least 5% relative to its 
previous size and that leads to a final disc size r disc <10 AU 
is denoted as ’’disc-destroying encounter 11 . Subsequent, poten¬ 
tially even stronger encounters are not considered in the anal¬ 
ysis because the disc is already destroyed. For these very small 
discs, effects such as viscous accretion and internal photoevap¬ 
oration may become important, but they were not considered 
in this work because of the computational effort. Therefore, all 
discs smaller than 10 AU are denoted ’’destroyed 11 . The term 
’’strongest encounter 11 in this work describes the encounter with 
the strongest influence on the disc, which is not necessarily the 
closest encounter, because the encounter mass ratio is also fac¬ 
tored in. 


Article number, page 4 of 


3. Results 


Previous studies of protoplanetary discs mostly investigated the 
effect of star-disc encounters on the disc mass (e.g. |Qlczak et al.| 
2010 Steinhausen & Pfalzner|2014|. Before studying this effect 


on the disc size in an entire star cluster, we first study the influ¬ 
ence of a fly-by on a single star-disc system. Recently, |Rosotti| 
|et al.[ ( |2014| > noted in their simulations of stellar clusters that 
more discs are affected by disc truncation than disc-mass los^j 
Here we wish to quantify this higher sensitivity of the disc size 
by systematically comparing it to disc-mass loss calculations. 
We used the results of Breslau et al. (2014) and Steinhausen et al. 
( |20 1 2[ > to obtain the disc sizes and masses after stellar encoun¬ 
ters. In contrast to the disc size, the disc mass-loss strongly de¬ 
pends on the distribution of material in the disc. Here, we used a 
disc-mass distribution that follows a power-law: S oc r 1 . 


3. 1. Sensitivity of disc properties 

In Fig.[T]a), the size (asterisks) and mass (open squares) of a 
disc after an encounter with an equal-mass perturber (mn = 1) 
are shown as a function of the relative periastron distance 
r p = r pe ri / r prev i ous , where r per i is the periastron distance and 
previous the disc size before the encounter. The effect of an en¬ 
counter on the size is much stronger than on the disc mass. For 
example, a close grazing encounter, equivalent to r p = 1, reduces 
the disc size to nearly a quarter of its previous size (28%), but it 
still retains more than half of its mass (57%). 

This higher sensitivity of the disc size compared to the disc 
mass also holds true for encounter partners that differ in mass 
(mi 2 + 1). The shaded areas in Fig. |T]b) depict the parameter 
space in which the encounters reduce the disc size (asterisks) and 
disc mass (open squares) by at least 5%. At any given mass ratio 
a reduction of the disc size by 5% is always achieved already 
with a much more distant encounter than is necessary for the 
same reduction in mass. And more importantly, a distant fly-by 


1 They used a disc-mass density distribution according to E oc r 2/3 
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distance to cluster centre [pc] 

Fig. 2: Median disc size as a function of the distance to the 
cluster centre for model D2 with different initial disc sizes r; n i t : 
100 AU (open triangles), 200 AU (filled circles), 500 AU (open 
circle), 800 AU (filled squares), 1 000 AU (open squares), and 
100000 AU (asterisks). 


can affect the disc radius without any change of the disc mas^] 
In other words, our data show that, for example, an encounter 
leading to 5% disc-mass loss can reduce the disc size down to 
60% of its previous value. This knowledge can now be used to 
quantify the sensitivity of the disc size compared to the disc mass 
in a whole cluster population. 


3.2. Initial disc size 

Most observed discs in clusters have sizes of 100 - 200 AU. 
However, a few protoplanetary discs with radii of a few 100 AU 
were found in the ONC as well. This indicates that a broad spec¬ 
trum of sizes exists. It is unclear whether these observed disc 
sizes are representative for their initial state or whether they 
changed during the 1-2 Myr since their formation, nor do we 
know which role the environment played in this context. As a 
consequence, there is no information available about a natal or 
initial disc-size distribution within the ONC, for instance. 

For this reason, we first investigated how the choice of an ini¬ 
tial disc size affects the final disc-size distributions in a cluster. 
We started our investigation with the model cluster D2, which is 
representative for the ONC. All stars in this cluster were initially 
surrounded by discs of the same size rj n j t . This initial disc size 
was then varied from 100 AU to 100 000 AU. Although most en¬ 
counters occur in the first Myr of cluster evolution ( |Pfalzner| 
et al .|2006 1 , the simulation time was chosen to be 5 Myr to en¬ 
sure that all relevant encounters are accounted for. 

To obtain an overview of the disc-size distribution in the 
cluster, the median disc sizes in different radial bins of the cluster 
were determined. Note that these median disc sizes were deter¬ 
mined for each distance-bin separately, not in a cumulative way. 
Figure [2] shows these median disc sizes for different initial disc 
sizes - 100 AU (open triangles), 200 AU (filled circles), 500 AU 
(open circle), 800 AU (filled squares), 1 000 AU (open squares), 
and 100000 AU (asterisks). 

Near the cluster centre the median disc sizes are significantly 
smaller than in the cluster outskirts. For example, for an initial 


2 The disc-mass change is also less sensitive than the angular momen¬ 

tum change, see|Qlczak et al.|f2006| l; |Pfalzner & Qlczak| ( [2007| > ; [Stein- 1 


hausen et al. (20121. 


Table 4: Fractions of discs influenced by encounters. 


Model 

Pcluster 

[ pc' 2 3 ] 

./reduced 

f< \ 000AU 

/<100AU 

/destroyed 

DO 

15 

65.1% 

42.8% 

8.8% 

3.4% 

D1 

30 

75.4% 

53.7% 

11.0% 

3.7% 

D2 

61 

82.7% 

64.9% 

14.6% 

4.3% 

D3 

122 

87.1% 

74.2% 

20.4% 

5.2% 

D4 

245 

89.2% 

81.3% 

28.8% 

6.4% 

D5 

489 

88.7% 

83.9% 

39.2% 

8.7% 


Notes. Fraction of stars with disc-size-reducing encounters (Col. 3), 
with discs smaller than 1 000 AU (Col. 4), with discs smaller than 
100 AU (Col. 5), and with destroyed discs (Col. 6) in the different clus¬ 
ter models (Col. 1) after 5 Myr of cluster evolution. The average number 
densities of the cluster models are given in Col. 2 (rounded values, taken 
from Table [3}. 


disc size of 500 AU, the median disc size at the rim of the cluster 
core (< 0.3 pc) is < 100 AU, but farther out at 2 pc it is more than 
five times as large (~ 500 AU). This trend might be expected be¬ 
cause the density is much higher in the inner parts of the clusters. 
Thus, more encounters occur there and more discs are reduced 
in size. Furthermore, the encounters in the inner cluster regions 
are on average closer, that is, they have lower periastrons (see 
e.g. Scally & Clarke||2001| >, which is due to the higher density. 
Therefore, the discs are smaller than in the less dense cluster 
outskirts. 


Overall, if the discs in the ONC were initially larger than 
100 AU they are significantly reduced in size by stellar encoun¬ 
ters, especially in the cluster core (< 0.3 pc). Any initial disc size 
Unit > 100 AU leads to the same median disc size distribution in 
the cluster core. In addition, an initial disc size > 1 000 AU was 
assumed; the final median disc size distribution in the cluster is 
the same up to a distance of ~ 4 pc, which in the case of the 
ONC would be the entire cluster (/fluster = 2.5 pc). This finding 
was then used to set up discs in the different cluster types to study 
their influence on the disc-size distribution. 


3.3. Numerical experiment 

Building on this finding that the median final disc size in a clus¬ 
ter is largely independent of the initial disc size, we performed a 
numerical experiment to be able to separate processed from un¬ 
processed discs. The initial disc sizes in each cluster were set 
t° Unit = 100000 AU (asterisks in Fig.|2]i. Although this is not 
a realistic scenario, it allowed us to determine the relevance of 
encounters in the cluster environments. 

This numerical experiment is applied to clusters of dif¬ 
ferent average densities, ranging from 15 pc~ 3 (model DO) to 
~ 500 pc~ 3 (model D5). For the densities of the other models, 
see Table [3 The cluster model D2 was denoted ONC model be¬ 
cause after 1 Myr of cluster evolution it was representative for 
today’s Orion nebula cluster, see Sect. 12.II Discs that were de¬ 
stroyed (rji sc < 10 AU) were excluded from the data set and 
treated separately. If not indicated otherwise, the outcome of the 
cluster simulations after 5 Myr of evolution is shown. 


3.3.1. Disc-size distributions 

As a first step, we quantified the amount of discs that are in¬ 
fluenced by encounters in the different cluster environments. 
Figure [3] shows the fraction of stars that undergo at least one 
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Fig. 3: Fraction of stars with destroyed discs (open circles) and 
whole discs (filled circles) as a function of the average cluster 
density. Depicted are discs smaller than 100000 AU (solid line), 
smaller than 1 000 AU (dashed line), and smaller than 100 AU 
(dashed-dotted line) after 5 Myr. 

disc-size reducing encounter (solid line, filled circles) as a func¬ 
tion of the average cluster number density. The fraction of stars 
with disc-size-reducing encounters increases, as expected, with 
cluster density. In the least dense clusters, some stars do not un¬ 
dergo any encounter at all or only very weak ones, that is, low 
mass ratios and high periastrons. Here the fraction of stars with 
disc-size-reducing encounters is ~ 65%. By contrast, in the dens¬ 
est cluster, stellar encounters are very frequent and stronger, that 
is, high mass ratios and/or lower periastrons. In this cluster type, 
nearly all stars (~ 90%) undergo at least one disc-size reducing 
encounter. 

The initial size of the discs was set to 100000 AU, thus, 
after the first disc-size-reducing encounter it might still be as 
large as 95 000 AU. But most encounters reduce the discs to 
much lower values and most discs undergo several disc-size- 
reducing encounters in the 5 Myr of cluster development. Con¬ 
sequently, two questions arise: i) what is the final disc-size distri¬ 
bution? and ii) how many discs are destroyed? The dashed and 
dashed-dotted lines with filled circles in Fig. [3] depict the frac¬ 
tions of discs with a final size smaller than 1 000 AU and 100 AU 
(see also Table [4]). The fraction of stars with discs smaller than 
1 000 AU nearly doubles from low-density clusters (43%) to 
high-density clusters (84%). Nevertheless, even in environments 
with densities of 15 pc~ 3 , stellar encounters significantly reduce 
the sizes of discs. The slope of the curve flattens towards the 
high-density cluster end, see Fig. [3] This flattening can be ex¬ 
plained by studying the fraction of destroyed d i sc ‘Fj( dotted line, 
open circles). It amounts to ~ 9% in model D5, which means 
that in this case, nearly all discs (97%) are either influenced by 
a disc-size-reducing encounter or are destroyed. This is near sat¬ 
uration, which is the reason for the flattening of the slope. Fur¬ 
thermore, the fraction of disc-destroying encounters triples from 
~ 3% in model DO to ~ 9% in model D5 (see Tabled, whereas 
the whole discs < 1 000 AU only double. This additionally limits 
the fraction of discs smaller than 1 000 AU. 

The fraction of stars with final disc sizes < 100 AU also in¬ 
creases with cluster density (dashed-dotted line). In the high- 
density cluster model it is four times higher than in the low- 

3 A detailed description of the properties of the disc-destroying en¬ 
counters can be found in the appendix. 
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distance to cluster centre [pc] 

Fig. 4: Median disc size in clusters of different densities (rising 
from DO, open squares, to D5, open triangles) for different dis¬ 
tances of the disc-hosting star to the cluster centre (bins). The 
line with asterisks indicates the ONC model, with the grey area 
containing 50% of all stars with disc sizes around the median 
disc size. For explanations, see text. 


density model, ~ 40% and ~ 10%. In this case, only 19% of 
all discs in model D5 are either reduced to sizes <100 AU 
or destroyed, so there is no flattening of the slope due to 
near saturation. The fraction of discs < 100 AU increases much 
more strongly with cluster density than the fraction of discs 
< 1 000 AU. For example, in the least dense cluster the ratio 
between discs smaller than 1 000 AU and discs smaller than 
100 AU is roughly 5:1, whereas for the densest cluster it in¬ 
creases to nearly 2:1. The reason is that the increasing cluster 
density leads to closer encounters that in turn reduce many more 
discs to smaller sizes. Summarising, the final disc size distribu¬ 
tion and the fraction of destroyed discs strongly depends on the 
average cluster number density. 

3.3.2. Distance to the cluster centre 

We have shown that the final disc sizes strongly depend on the 
average number density of their cluster environment. Each of 
these clusters was set up with a density profile according to 
EqCl so they are much denser in the inner regions than in the 
outskirts. Therefore, the disc sizes should also be a function of 
the position of the disc-hosting star in the cluster. In Fig. [4] we 
compare the median disc sizes as a function of the star’s final 
(after 5 Myr) distance to the cluster centre in the different clus¬ 
ters. For comparison, the line with asterisks is the same as in 

Fig-0 

As expected, the median final disc size increases with the dis¬ 
tance to the cluster centre for all average cluster densities. Stars 
that are finally located at a distance of 0.3 pc from the cluster 
centre have a median disc size of about 200 AU in the sparsest 
cluster (DO, open squares), whereas in the cluster outskirts it is 
much larger, for instance ~ 2 400 AU at a distance of 3 pc. For 
the densest clusters model D5, the median disc sizes are much 
smaller as a result of the elevated rate of close encounters, as 
described before, but the dependence on the host star position in 
the cluster is the same. At 0.3 pc the median disc size is about 
20 AU, whereas at 3 pc it is ~ 300 AU (open triangles). 

To obtain the variation of disc sizes in the clusters, the shaded 
area in Fig. [4] depicts the disc-size range for 50% of all stars in 
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spectral type 
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Fig. 5: Median disc size in cluster models DO (open squares), D2 
(asterisks), and D5 (open triangles) as a function of the stellar 
mass [ M q \ and spectral type. 


the ONC model that have disc sizes of about the median size 
(asterisks). For example, half of the stars with final distances of 
0.3 pc to the cluster centre, which represents the border of the 
Trapezium in the ONC, have discs with sizes of 40 - 200 AU. 
The median disc size at this distance is 100 AU. At the cluster 
radius of the ONC (2.5 pc), half of all discs have sizes between 
300 - 1 300 AU, the median disc size is 630 AU. This means that 
half of the stars in the ONC model have disc sizes that vary by 
about a factor of two about the median, but which are still much 
smaller than the initial disc size of 100000 AU. The final disc 
size is a function of the host star position in the cluster, or in 
other words, not only of the global average cluster density, but 
also of the local density of its environment. 


3.3.3. Disc sizes for different stellar types 

As mentioned in Sect.Q] it is still a subject of debate whether 
the disc size is a function of the stellar mass. We investi¬ 
gated whether encounter-induced disc-size reduction might be 
responsible for this potential dependency. Figure [5] shows the 
median disc size as a function of stellar mass and spectral 
type. For clarity, only the results for models DO (open squares), 
D2 (asterisks), and D5 (open triangles) are shown. 

The median disc size does indeed depend on the stellar mass. 
For M- to B-type stars in the ONC model, the median disc size 
increases by roughly a factor of two, from 420 AU to 880 AU. 
Around O-stars the discs are smaller than around all other star 
types with a median size of 270 AU. The same trend in me¬ 
dian disc sizes can be found for models of lower density (e.g. 
model DO in Fig. [5j. The smaller sizes of discs around O-stars 
in models DO and D2 can be explained by the encounter fre¬ 
quency, which as a function of the relative perturber mass has 
two peaks (see Fig. 4 of Olczak et al.|2010) >. One peak lies at low 
mass-ratios, the other one at high mass-ratios. For model DO the 
high-mass stars dominate the encounter history because of grav¬ 
itational focusing. Therefore, the difference between the median 
disc sizes of M-stars and O-stars is largest in this model. 

In contrast, in model D5 the median disc size of O-stars 
(130 AU) is slightly larger than the one of M-stars (~ 90 AU). In 
this cluster type, the low-mass stars dominate the encounter his¬ 
tory so that M-stars have slightly smaller discs than all other star 
types - including O-stars. The ONC model D2 represents a turn¬ 


ing point: the most massive and least massive stars are equally 
important as encounter partners. Thus, the absolute difference 
between the discs around M-stars and O-stars is smaller than in 
model DO. 

All cluster models were set up without primordial mass seg¬ 
regation. Nevertheless, in the 5 Myr of cluster evolution, dynam¬ 
ical mass segregation might occur. The final positions of the 
B- and O-stars showed no mass segregation in the low-density 
clusters and only weak mass segregation in the densest models. 
Therefore, the impact of gravitational focusing on the disc size - 
especially in the low-density clusters - is most probably stronger 
than the effect of dynamical mass segregation. 


4. Discussion and conclusion 


Our simulations include several assumptions. Their influence on 
the results is discussed below. 

Most stellar clusters show signs of mass segregation, but 
whether this is the case for the ONC is still under debate, cf. 
|Bonnell & Davies] ( |1 998) i; |Qlczak et al.| ( |201 1) ; |Huff & St abler 
(2006). It is currently unclear whether clusters are primordially 
mass-segregated or not. Here we did not take into account any 
primordial mass segregation, only dynamical mass segregation 
during the simulation. This was done for simplicity because we 
aimed to study clusters of different densities and their influ¬ 
ence on disc sizes. Primordial mass segregation would have in¬ 
creased the initial cluster-core density and altered the encounter 
frequency. Primordial mass segregation affects this type of clus¬ 
ters in two ways: i) the most massive stars undergo many more 
encounters as a result of their location in the dense cluster centre, 
and therefore the disc sizes are smaller around O-stars than in the 
case without mass-segregation; ii) the stronger gravitational fo¬ 
cusing increases the number of encounters. This leads to a some¬ 
what lower disc frequency and smaller disc sizes, mainly around 
the O-stars. 

All stars in the clusters were assumed to be initially sin¬ 
gle, neglecting primordial binaries. If these were included, more 
discs would be truncated to a smaller size or even destroyed by 


the binary stars, especially for very close binaries. Kohler et al. 
( |2006| found that the binary frequency in the ONC is higher for 
stars with masses > 2 M 0 than for low-mass stars. Generally, 


the binary frequency increases with stellar mass, see Duchene 
|& Krausj ( 2013| ). In our simulations, primordial binaries would 
therefore lead to even smaller discs for B- and O-stars than an¬ 
ticipated in Fig. [5] 

The cluster models investigated are representative for young, 
massive clusters that are still embedded in their natal gas. In the 
solar neighbourhood the star formation efficiency is relatively 
low (30%) and the gas expulsion drives the expansion of the 
cluster, while at the same time a large portion of its stars be¬ 
comes unbound ( |Pfalzner & Kaczmarek|2013) . Observations of 
stellar clusters indicate that after 1-3 Myr the gas is expelled 
from the cluster and thus the star-formation process ends. We 
did not consider the effect of gas expulsion. In a follow-up pa¬ 
per, this case of cluster dynamics will be studied in detail. 

We chose an artificially large initial disc of 100000 AU 
around every star to determine the relevance of encounters in a 
cluster environment. Such discs are not physical, but have to be 
seen as part of the performed numerical experiment. The result¬ 
ing median disc sizes can thus not represent any real distribution 
as found in the ONC, for example, but rather give a first estimate 
of the influence of encounters on discs after 5 Myr of evolution, 
especially when considering that other effects are neglected (see 
below). The surprisingly good match to the observed data shows 


Article number, page 7 of 10 































A&A proofs: manuscript no. strong_effect_of_cluster_environment_on_protoplanetary_disc_sizes 


that the process of disc-size reduction that is due to stellar en¬ 
counters is indeed important and should not be neglected. The 
largest discs in the ONC could have formed even larger and al¬ 
ready have been pruned to their current size by dynamical inter¬ 
actions. 

All encounters in the clusters were assumed to be prograde, 
coplanar, and parabolic (eccentricity = 1). These are the most 
destructive encounters (Clarke & Pringle 1993] Pfalzner et al. 


2005b). Especially in very dense clusters, a large portion of en¬ 
counters can have eccentricities > 1. Such eccentric orbits are 
less efficient in removing disc mass (Olczak et al. (2010|>). In ad 


dition, inclined encounters are common in cluster environments. 
Coplanar encounters lead to smaller disc sizes than their inclined 
counterparts. Therefore, the cases of parabolic coplanar encoun¬ 
ters should be regarded as lower limits (Breslau, in prep.). 

There are other effects that might influence the disc size. 
First, we used the formula of [Breslau et al.| ( j2014| ) to de¬ 
scribe star-disc encounters. It might have been more appropri¬ 
ate to investigate disc-disc encounters. However, it was shown 
by|Pfalzner et al.| (2005a)> that in disc-disc encounters between 


low-mass discs (0.01 M Q ) material is exchanged, but it is mainly 
transferred to the inner disc regions. Since the disc size is de¬ 
termined using the steepest point in the outer slope of the disc 
density distribution, this captured disc material is not expected 
to alter the final size much. 

Second, it has been suggested that photoevaporation might 
be an important mechanism to shrink or completely destroy discs 


(see e.g. Storzer & 

iollenbach 1999\ 

Scally & Clarke 

2002 

lohnstone et all2004 

Adams et al. 2006-Alexander et al. 

2006 


2009) ). However, we modelled stellar clusters in the embedded 


phase, during which stage the presence of the natal gas is ex¬ 
pected to reduce the effect of photoevaporation. If photoevapo¬ 
ration nevertheless played a role, the presented disc sizes would 
be even smaller for stars close to the most massive stars, espe¬ 
cially near the cluster core. 

Third, the discs were modelled consisting of mass-less tracer 
particles and hydrodynamical effects were excluded, for an in- 
depth discussion see |Breslau et al,|f2014[ ). |Rosotti et al 7| ( |20 1 4[ ) 
found that the hydrodynamical evolution of the discs is impor¬ 
tant for high viscosity. In this case, the discs spread quite rapidly 
after an encounter and compensate for the size reduction that is 
due to encounters. In the case of low viscosity, it was the closest 
encounter alone that determined the final disc size. The encoun¬ 
ters were particularly important in the dense central cluster areas. 
We did not include viscous disc spreading in our simulations. If 
we had considered this effect, the discs that already underwent 
their strongest encounter (high mass ratio and/or small perias- 
tron distance) would have been most affected by changes due to 
disc spreading. If the disc had had time to expand again, the fi¬ 
nal disc size after a subsequent relatively strong encounter would 
have been larger than presented here. To study this effect in more 
detail, the next natural step would be to include viscous spread¬ 
ing in our type of simulations. 

As a result of its sensitivity, the disc size is a suitable tool 
for quantifying the influence of the stellar environment. One in¬ 
dication that this might be the case comes from the observations 
by Vicente & Alves ( |2005| ). They found that roughly 60% of the 
Trapezium stars (< 0.3 pc) have discs smaller than 50 AU. Our 
simulations indicate that for model D2 (ONC), the median disc 
size of all stars in the cluster core (< 0.3 pc) is roughly 80 AU. 
This fits the observed values remarkably well, especially when 
considering that our simulations only take into account stellar 
encounters as a source of disc-size influence. They play an im¬ 


portant role in terms of disc-sizes but should be taken together 
with the other effects discussed above (photoevaporation, vis¬ 
cous disc-spreading, etc.) to reproduce the actual disc-size dis¬ 
tributions as found in stellar clusters today. 

The protoplanetary discs may eventually form planetary sys¬ 
tems. Their potential physical extent, as shown above, is defined 
to a high degree by the birth environment. The influence of a stel¬ 
lar cluster on already formed planetary systems has been studied 
theoretically in the past (e.g. Craig & Krumholz 2013; Malmberg 
|et al.|201~T| ). |Craig & Krumholz| ( |2013| ) found that in dispersing, 
substructured clusters planets close to their host stars (few tens 
of AU distance) are not perturbed by stellar flybys. Compared to 
the spherical clusters modelled here, the substructures lead to re¬ 
gions of enhanced stellar density and thus to a higher encounter 
frequency in certain areas of the cluster. These substructures dis¬ 
solve in most cases at the latest when the whole cluster disperses. 
When the substructures are taken as regions of (transiently) en¬ 
hanced stellar density, our results support their findings, as only a 
few per cent of discs are completely destroyed (< 10 AU), while 
most discs are finally larger than a few tens of AU, see Fig. [3] 

5. Summary 

We investigated the influence of the cluster environment on the 
disc size. Combining Nbody simulations of clusters of different 
densities with simulations of protoplanetary discs after star-disc 
encounters, we found that unlike the disc mass, the disc size 
can already be changed by relatively distant encounters. We 
performed a numerical experiment, choosing a very large initial 
disc size around each star. The results show that encounters 
are capable of reducing disc sizes of most stars to values 
< 1 000 AU in all investigated clusters. The resulting sizes agree 
well with typical disc sizes found by observations in the ONC 
today, for example. It might be that the cluster environment 
rather than the initial size determines the final extent of a 
protoplanetary disc in many cases. 


Appendix A: Properties of disc-destroying 
encounters 

Figure [AT] a) shows the fraction of disc-destroying encoun¬ 
ters as a function of the encounter mass ratio for three 
cluster models (DO, D2, D5). For the densest cluster model 
(D5, open triangles), the encounter mass-ratio distribution has 
two peaks: one at low mass ratios (mn = 3) and a second smaller 
peak at very high mass ratios {mn = 250). So the disc-destroying 
encounters most probably have quite low mass ratios in this en¬ 
vironment. This is in accordance with the results of lOlczak et al.l 
( |2010] >. They found that in clusters of higher density than the 
ONC the low-mass stars dominate the encounter history, lead¬ 
ing to lower encounter mass ratios. They stated that the ONC 
itself is a turning point, low-mass and high-mass stars similarly 
contribute to the encounter history. This is reflected in the en¬ 
counter mass-ratio distribution for the ONC model (asterisks), 
which is flat and does not show any sign of a preferred mass 
ratio for disc-destroying encounters. For the least dense cluster 
model (DO, open squares), one would expect, based on the work 
of |Qlczak et al. ( 2010| ), that high mass-ratio encounters should 
be largely responsible for disc destruction, as the most massive 
stars undergo many encounters due to gravitational focusing. But 


this is not the case, as Fig. A.l a) shows. The distribution of 
disc-destroying encounters has a very broad peak at low mass 
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Fig. A. 1: a) Fraction of disc-destroying encounters as function of the mass ratio and b) fraction of disc-destroying encounters as 
function of the periastron distance for cluster models DO (open squares), D2 (asterisks), and D5 (open triangles). 


ratios. This might be due to the mass sampling in this work: the 
masses were randomly sampled from an IMF (see Sect. [2J so it 
is more probable to obtain very massive stars in clusters with 
higher membership. In low-mass clusters the most massive star 
is - on average - less massive than in the high-mass cluster. 

Figure |AJ~| b) shows the fraction of disc-destroying encoun¬ 
ters as a function of their periastron distance for the same clus¬ 
ter models as above (DO, D2, D5). For all cluster types most 
disc-destroying encounters typically have distances of about 
30 - 40 AU. The distribution for model D2 is slightly broader 
than those of the other cluster models. The disc size is deter¬ 
mined by both the periastron distance and the mass ratio (see 
Eq-0, thus, this is in accordance with the broad spectmm of 
mass ratios for the disc-destroying encounters in Fig. |A.l| a). If 
the mass ratio is low, the encounter has to be closer, meaning 
that the periastron distance needs to be smaller, to destroy a disc 
than is necessary in a high mass-ratio case. 


Appendix B: Location of disc-destroying 
encounters 

The larger part of discs is influenced by their stellar encounters 
(65 - 89%, see Table [4}. but is not destroyed. As already men¬ 
tioned in Sect. 13.3.11 in the sparsest cluster, about 3% of the discs 
are destroyed after 5 Myr; for the densest cluster this fraction in¬ 
creases to ~ 9%. This means that the fraction of destroyed discs 
increases with increasing average cluster density. Analogously 
to the analysis in Sect. 13.3.21 we can now determine whether the 
disc-destroying encounters are also a function of the star’s po¬ 
sition in the cluster, or in other words, a function of the local 
density. 

Figure |BTT| shows the position of the disc-hosting stars at 
the time of the encounter for three different cluster models 
(DO, open squares, D2, asterisks, and D5, open triangles). For 
the least dense cluster, most disc-destroying encounters take 
place at a distance of 0.1 - 0.2 pc from the cluster centre. In 
the case of the ONC model cluster (D2, asterisks) the whole 
curve is shifted towards lower distances. As the cluster is much 
denser than model DO, more destructive encounters occur at such 
small distances from the cluster centre. For the densest cluster 
model (D5, open triangles), the distribution is entirely different. 
In this cluster type, the encounters are dominated by low-mass 
stars, in contrast to the low-density clusters where the most mas- 



distance to cluster centre [pc] 

Fig. B.l: Fraction of disc-destroying encounters as function of 
distance to the cluster centre where they took place for cluster 
models DO (open squares), D2 (asterisks), and D5 (open trian¬ 
gles). 


sive stars dominate the encounter history ( Olczak et al.||2010| >. 
As a result, the distribution of locations is much shallower and 
peaks only in the most inner part of the cluster (~ 0.01 pc). 
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